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Even after 25 years of research the pairing mechanism and - at least for electron doped com- 
pounds - also the order parameter symmetry of the high transition temperature (high-Tc) cuprate 
superconductors is still under debate. One of the reasons is the complex crystal structure of most 
of these materials. An exception are the infinite layer (IL) compounds consisting essentially of 
Cu02 planes. Unfortunately, these materials are difficult to grow and, thus, there are only few 
experimental investigations. Recently, we succeeded in depositing high quality films of the electron 
doped IL compound Sri_xLaxCu02 (SLCO), with x ~ 0.15, and on the fabrication of well-defined 
grain boundary Josephson junctions (GBJs) based on such SLCO films. Here we report on a phase 
sensitive study of the superconducting order parameter based on GBJ SQUIDs from a SLCO film 
grown on a tetracrystal substrate. Our results show that also the parent structure of the high-Tc 
cuprates has (i3,2_^2-wave symmetry, which thus seems to be inherent to cuprate superconductivity. 



PACS numbers: 74.50.+r, 74.72.-h, 74.72.Ek, 85.25.Dq 



Since the discovery of high transition temperature 
(high-Tc) superconductivity in cuprates^, tremendous 
work has been performed on these materials. Researchers 
succeeded in increasing Tc from initially 30 K to 135 K^^ 
by synthesizing increasingly complex compounds. How- 
ever, the microscopic mechanism causing high-Tc super- 
conductivity still has not been identified and is one of 
the biggest issues in solid state physics. All these ma- 
terials have in common that superconductivity resides 
in the copper oxide (CUO2) planes where superconduct- 
ing charge carriers, Cooper pairs, form. An "infinite 
layer" (IL) cuprate consisting essentially of Cu02 planes 
is therefore of fundamental interest for all questions ad- 
dressing the basics of high-Tc superconductivity. In 1988, 
Siegrist et al. succeeded in synthesizing such a simple 
cuprate, which is known as the "parent structure" of 
cuprate superconductors^. Its Cu02 planes are only sep- 
arated by a single alkaline earth metal plane {A = Ca, 
Sr or Ba), forming a ACUO2 crystal. Upon electron- 
doping, it turned out to be superconducting with maxi- 
mum Tc= 43 K^^ 

A striking and highly debated feature of the cuprate 
superconductors is their unconventional order parame- 
ter symmetry. Whereas for hole-doped cuprates dx2_y2- 
wave pairing has been established^^^^, for electron- 
doped cuprates ^^ the issue is not yet completely set- 
tled. After controversial discussion, the electron-doped 
T'-compounds T2-a^Cea^Cu04 (T = La, Pr, Nd, Eu 
or Sm) have been shown to be predominant dj.2_y2- 
wave superconductors by a number of phase-sensitive 
experiments^^^^. By contrast, for the parent com- 
pounds the pairing symmetry is essentially still unknown, 
since a variety of experimental tests yielded conflicting 
results^^'^^^^. Phase-sensitive tests, such as experiments 
on corner junctions^^, tricrystal rings^^ or tetracrys- 



tal SQUIDs^^'^^'^^, are widely recognized to provide a 
clear evidence for the pairing symmetry of the order 
parameter ^^. Such experiments rely on Josephson junc- 
tions, which for IL cuprate thin films became available 
only very recently^^. 

Here we report on the fabrication and characterization 
of thin-film Sri_a^Laa;Cu02 (SLCO) dc SQUIDs based 
on BaTiOs-buffered tetracrystal SrTiOs substrates. The 
geometry involved is designed to be frustrated for 
(i3,2_^2-wave pairing, i.e. the SQUID ring comprising the 
tetracrystal point contains one junction and one tt junc- 
tion, if the order parameter has (ia,2_^2-wave symmetry. 
This device will be referred to as 7r-design SQUID. Its ge- 
ometry, together with the design of a reference SQUID, 
is sketched in Fig. 1. There are four grain boundaries 
(GBs), labelled 1-4. GBs 1-3 have misorientation an- 
gles of 30°, while the misorientation angle of GB 4 is 
0°. The TT-design SQUID comprises all GBs. GB 4 will 
not form a grain boundary junction (GBJ) due to its 0° 
misalignment angle, in contrast to GBs 1-3. GBs 2 and 
3, having a width of 58 /im, form the active Josephson 
junctions in the current and voltage lead configuration 
indicated in Fig. 1. The bias current also passes GBJ 1 
which, however, is much longer (~1.5mm) than GBJs 2 
and 3 and thus will have a much higher critical current. 
Below, we will see however, that flux quanta (Joseph- 
son fluxons) can enter this GBJ, which thus cannot be 
ignored in the data analysis. If SLCO is a (ia,2_^2-wave 
superconductor, one of GBJs 2 and 3 faces a sign change 
of the order parameter (GBJ 3 in Fig. 1), thus forming a 
TT Josephson junction. The other GBJs are conventional. 
The area of GBJs 2 and 3 is not much smaller than the 
area of the SQUID hole. In this "spatially distributed 
junction" design^^ the junction's Ic vs. H modulation 
(Fraunhofer pattern) is superposed on the SQUID mod- 
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FIG. 1: (Color online) Schematic layout of the SQUIDs. The 
0-design SQUID comprises two conventional GBJs (0 junc- 
tions) straddling a single 30° [001]-tilt grain boundary. The 
TT-design SQUID comprises four GBs. The misorientation an- 
gle of GB 4 is 0°. All other misorientation angles are 30°. 
The dx2_y2-wave order parameter is indicated by the clover- 
leaf structure consisting of white and black lobes, indicating 
the sign change of the order parameter. Leads for bias cur- 
rent / and voltage V are indicated. In some experiments, for 
the TT-design SQUID we have also sent a current /mod across 
GB 1. Magnetic fields have been applied perpendicular to the 
substrate plane. 



ulation on a similar field scale. The (a)symmetry of the 
SQUID modulation relative to the Fraunhofer envelope 
allows to detect residual fields and often also trapped 
magnetic flux. 

The reference SQUIDs - there were two reference 
SQUIDs, producing very similar results - cross only one 
of the 30° GBs and incorporate two 50 fim. wide GBJs, 
which act as conventional junctions both for s-wave and 
dx2_y2-wdive order parameters. Below, these devices will 
be referred to as the 0-design SQUIDs. Both the TT- 
design SQUID and the reference SQUIDs had rectangular 
SQUID holes with an area As = 50 x 75 /im^. 

The samples have been fabricated by pulsed laser de- 
position, as described elsewhere^^'^^'^^. In brief, we first 
deposited a 25 nm thick BaTiOs thin-film on the SrTiOa 
tetracystal, acting as a buffer layer. This layer was fol- 
lowed by a 22 nm thick SLCO thin-film, with doping 
X ^ 0.15. Finally, a 10 nm thick gold layer was evap- 
orated in-situ, protecting SLCO from degradation and 
acting as resistive shunt for the GBJs. The SQUIDs 
were patterned by standard photolithography and argon 



ion milling. The SLCO film had a critical temperature 
Tc ^ 18 K. Electric transport measurements were per- 
formed at T = 4.2 K in a four-point configuration, with 
the sample mounted inside a noise-filtered, magnetically 
and radio frequency shielded probe in a liquid-helium de- 
war. A SQUID amplifier was used to allow for low- noise 
measurements. 

Below we discuss data of one of the 0-design SQUIDs 
and of the TT-design SQUID. The current voltage {IV) 
characteristics of these devices were nonhysteretic and 
could be well reproduced by the SQUID Langevin equa- 
tions, extended by taking the nonzero junction width into 
account. Details can be found in Ref. 36, giving rea- 
sonable values for the junction parameters Iq (maximum 
amplitude of Josephson current), R (junction resistance) 
and C (junction capacitance). The junctions of each de- 
vice were symmetric in terms of /q? R and C, with values 
/o = 8.2 /iA, R = 0.92 (^, C = 24pF (0-design SQUID) 
and /o = 12.2 /iA, R = 0.871], C = 25pF (TT-design 
SQUID). For the inductance parameter Pl = 2IqL/^q^ 
where <l>o is the magnetic flux quantum and L is the total 
inductance of the SQUID, we found ^l = 1-4 (0-design 
SQUID; L = 177 pH) and 2.2 (TT-design SQUID; L = 
187 pH), with an asymmetry a^ = 0.05 between the left 
and right arm of the SQUID (both designs). A fraction 
fj = 0.128 (0-design SQUID) and fj = 0.12 (TT-design 
SQUID) of the fiux ^ applied to the SQUID loop was 
coupled to each junction. The numbers for fj were de- 
rived from an analysis of the SQUID critical current Ic 
vs. applied field H (see below). 

Fig. 2 (a) shows the measured Ic vs. H for the 0- 
design SQUID (black line). Ic was determined using a 
voltage criterion Vc = 50 nV. To identify magnetic hys- 
teresis effects, /c vs. H was traced from 2.8 /iT to -2.8 /iT 
and back to 2.8 /iT. One observes a SQUID modulation 
with period jiq/S^H = 0.175 /iT on top of a Fraunhofer-like 
modulation, which is due to the finite junction size. No 
hysteresis is visible. The modulation period corresponds 
to an effective SQUID area of 1.18 • 10^ /im^, pointing 
to a flux compression of about 3.15. This is reasonable 
for our structures^^'^^. The insets of Fig. 2 (a) show Ic 
vs. H near H = {) ioi both positive and negative I^ 
The Ic maximum is close to ilf = 0, with a small off- 
set of 10.6 nT for positive Ic and 5.2 nT for negative I^ 
We contribute the asymmetry of ±2.7 nT in offsets to an 
inductance asymmetry (a/, = 0.05) of the two SQUID 
arms and the average part of 7.9 nT to residual fields in 
the cryostat. The red line in Fig. 2 (a) is a numerical cal- 
culation. It produces data very well inside the main max- 
imum of the Fraunhofer envelope. The first Fraunhofer 
side-maximum is lower in amplitude than the experimen- 
tal data, presumably due to a field distribution inside the 
junction, which is more complex than the homogeneous 
flux density assumed in our model. Most importantly, 
however, we see that the 0-design SQUID behaves as it 
should be expected from a conventional 0- SQUID. 

The measured /c vs. H of the TT-design SQUID is 
shown by the black line in Fig. 2 (b). Also here we have 
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FIG. 2: (Color online) Critical current Ic vs. applied mag- 
netic field H of (a) the 0-design SQUID and (b) the vr-design 
SQUID (black lines) together with calculated curves (red 
lines) for (a) a SQUID and (b) a tv SQUID. Insets show 
Ic vs. H on expanded scales. A 50 nV voltage criterion was 
used to determine Ic. 



varied H from 2.8 /iT to -2.8 /iT and back to 2.8 /iT. Ic 
is at a minimum near i7 = - a feature which appears 
when one of the two Josephson junctions exhibits an ad- 
ditional 7T shift in its phase. At negative Ic the minimum 
is at jiioH ^ 0.15 nT, while at positive Ic it appears at 
/JjqH ^ 7.8 nT, pointing to an offset field of about 4nT 
and a small asymmetry in inductance (a^ = 0.05). The 
SQUID modulation period is fioAH = 0.136 /iT, corre- 
sponding to an effective area of 1.52 • 10^ /im^ and a flux 
compression factor of 4.05. The overall modulation of Ic 
vs. H is described reasonably well by numerical calcula- 
tions (red line), however less well than /c vs. H of the 
0-design SQUID. 

A prominent feature are the jumps in Ic, visible at 
jiioH > 89 nT at positive Ic and at jj^qH < -80 nT at neg- 
ative Ic. There is only a very tiny hysteresis associated 
with these jumps, which is not even visible in Fig. 2 (b). 
By comparing measured and calculated Ic vs. H curves 
within the main Fraunhofer lobe one sees that the calcu- 
lated curve exhibits one additional SQUID period. These 
features indicate that magnetic flux quanta enter the de- 
vice at each jump. This effect was not visible for the 0- 
design SQUID. A strong candidate for flux entry is thus 
GBJ 1 which is absent in the reference SQUID. Note that 
the Ic jumps visible in Fig. 2 (b) occur point- symmetric, 
i.e., at positive Ic they occur at positive fields, while at 
negative Ic the field is negative, with about the same am- 
plitude as for positive Ic. This feature can also clearly 
be seen in ^ vs. <l> patterns taken at many values of bias 



22 


.(a) 


A A 


T 


= 4.8K 
= 2fjV 


20 


A 


/ A 


fei 


- 


_" 18 


(1) / \ 


1 y 


r 


y- 


16 


\ / 


V M / 


f\ 


\j 




V 


(2) (3) (4) 


(6) 


w 


14 


1 




. 


IV _ 



-0.6 



-0.4 -0.2 

MqH (mT) 



0.0 



(b) 

(1) (2) (3) (4) (5) (6) (7) 




FIG. 3: (Color onUne) LTSEM data: (a) Ic vs. H, as mea- 
sured in the LTSEM setup (black and green lines distinguish 
sweep directions), (b) SV images, taken along GBJ 1 at / = 
17 /iA [horizontal line in (a)] at the magnetic field values (1) 
- (7) indicated in (a). Graph (c) shows a 6V image of the 
SQUID hole and GBJs 2 and 3 at / = 22 /iA and /j^qH = 
-0.19 /iT. This image has been superposed to scale to image 
(b)-(2) in order to indicate the position of the SQUID hole. 
For each image ^Umax is indicated inside the graphs. 



current, see Ref. 36. There we also show, that the point 
symmetry in /c vs. H holds even for large values of H^ 
and that applying an additional current Imod across GBJ 
1 alters the values of ilf in a way that is compatible with 
the notion of Josephson fluxons having entered GBJ 1. 

For a final proof we have imaged the current distri- 
bution of the TT-design SQUID using low temperature 
scanning electron microscopy (LTSEM). Details of the 
method can be found in Refs. 17,38. In brief, the pulsed 
electron beam, which is scanned across the sample, causes 
local heating by ~ IK. The measured integral quantity is 
the voltage V across the SQUID, which is biased slightly 
above Ic. The electron beam causes a change SV{x,y) 
depending on the beam position {x,y). When the beam 
is scanned across GBJs 2 and 3 near a Ic maximum, a 



positive signal appears, because Ic is lowered, causing 
a slight increase of V (cf. Fig. 3 (c)). When GBJ 1 
is free of vortices no signal is expected from this GB. 
By contrast, when Josephson fluxons are present, local 
heating will alter the screening currents around the flux- 
ons. In a heated area the Cooper pair density and thus 
the maximum supercurrent density is suppressed, caus- 
ing an increase of the Josephson length Xj and, thus, 
the fiuxon is virtually deformed towards the heated area. 
When the electron beam is between the fiuxon center and 
the SQUID hole this causes an increase of the fluxon's 
stray flux coupled to the SQUID and thus a change in 
Ic. In the opposite case the stray flux is decreased. A 
fiuxon will thus appear as a bipolar signal SV{x^ y)^ with 
increased/decreased voltage relative to the unperturbed 
value of V . In a similar way, Abrikosov vortices trapped 
in a YBa2Cu3 07 SQUID^^, as well as Josephson fluxons 
having entered a GBJ^^, have been imaged by LTSEM. 
Fig. 3 (a) shows Ic vs. H^ as measured in the LTSEM 
setup at T = 4.8 K using a voltage criterion Vc = 2 /iV. 
There is a stronger offset field (~ 0.25 /iT) than in the 
transport setup, causing a shift by about 2 SQUID modu- 
lation periods. Due to the Fraunhofer envelope this shift 
is straightforward to recognize. The two sweep directions 
of H have been distinguished by black and green lines. 
There is a nonhysteretic region around the offset field; at 
larger values of H^ Ic jumps occur, leading to magnetic 
hysteresis. SV images of GBJ 1, cf. Fig. 3 (b), have been 
taken at / = 17 jaA at the field values indicated by la- 
bels (l)-(7) in Fig. 3 (a). No contrast appears when the 
sample is biased near the offset field [images (2), (3)] or 
at a field smoothly extending this /c vs. H pattern [im- 
ages (4), (5)]. By contrast, when Ic jumps have occurred, 
we observed a periodically modulated signal, having a 
period decreasing with the field amplitude H relative to 
the offset field [images (1),(6),(7)]. This is very indicative 
of Josephson fluxons having entered GBJ 1. In Fig. 3 (c) 
we also show a SV image of the area around the SQUID 
hole and GBJs 2 and 3. The image has been taken in 
a separate run, because we did not want to disturb the 
images of Fig. 3 (b) by scanning across these GBJs. It 
has been taken at the maximum of the SQUID modu- 
lation at / = 22 /iA and fioH = -0.19 /iT. The image is 
superposed to scale to image Fig. 3 (b)-(2), to give an 
impression of the position of the SQUID hole and GBJs 
2 and 3 relative to the images of Fig. 3 (b). 



The LTSEM data clearly show that the /c vs. H region 
of interest is free of trapped flux; we thus feel safe in 
interpreting the 7r-design data in favor of a dx2_y2-wdive 
symmetry of Sri_a^Laa^Cu02. 

One may in addition ask about subdominant order pa- 
rameters. A real superposition dx2_y2±s is not very 
likely due to the tetragonal crystal symmetry but, if 
present, could lead to an asymmetry of the critical cur- 
rents of GBJs 2 and 3 and, in consequence, to a similar 
shift as the one we interpreted in terms of an inductance 
asymmetry. On the other hand, the 0-design SQUID 
should not show this asymmetry and we thus believe that 
an inductance asymmetry is more likely. By contrast, a 
complex admixture of a subdominant order parameter 
would lead to a ground state phase different from tt or 
0. Then, the SQUID modulation would shift relatively 
to the Fraunhofer envelope, making the amplitude of the 
inner Ic maxima asymmetric. This effect is not observed 
at least on a ~ 5% level. 

In summary, our data clearly show, that the supercon- 
ducting order parameter of the electron doped infinite- 
layer high-Tc cuprate Sri_a;LaajCu02 has (ia,2_^2-wave 
symmetry. The phase sensitive configuration used was 
a TT SQUID patterned on a tetracrystal. The parasitic ef- 
fect of Josephson fluxons entering one of the grain bound- 
ary junctions has been ruled out by direct imaging of the 
local supercurrent contribution. Sri_a^Laa;Cu02 has the 
most simple crystal structure of all high-Tc cuprates. We 
conclude that the dx2_y2-wdive symmetry is inherent to 
cuprate superconductivity and neither restricted to hole 
doping nor related to the complex crystal structures that 
complicates an analysis of almost all other cuprate su- 
perconductors. 
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In this supplement we provide additional experimental 
data for the 0-design SQUID and the 7r-design SQUID: 
(I) current voltage {IV) characteristics, (II) voltage V 
vs. applied field H^ (III) symmetry considerations on 
/c vs. H^ and (IV) a measurement of the current step 
height caused by the SQUID LC resonances. Most mea- 
surements are accompanied by simulations. The model 
equations are described in section I. 



I. CURRENT VOLTAGE CHARACTERISTICS 
AND MODEL 

We start to characterize our devices by discussing their 
current voltage {IV) characteristics. Fig. SI (a) shows 
two IV characteristics of the 0-design SQUID. The mag- 
netic field H^ applied perpendicularly to the substrate 
plane, was adjusted such that the SQUID's critical cur- 
rent was at its first (counted from zero applied field) 
maximum (black circles, jiqH = 0.01 /iT) or at its first 
minimum (grey circles, iiqH = 0.1 /iT). In both cases 
the IV characteristics were nonhysteretic. They can be 
very well fitted by numerically calculated curves, solving 
the SQUID Langevin equations^, extended by taking the 
nonzero junction width into account. The model assumes 
that the Josephson junctions can be described by the re- 
sistively and capacitively shunted junction model^'^. We 
have assumed further, that a fraction fj of the flux <l> 
through the SQUID loop is homogeneously coupled into 
the junctions, causing a linear increase of the Josephson 
phase differences 7/c {k = 1,2) inside the junctions. This 
contribution has been integrated out analytically, yield- 
ing a sine function for the critical current /q vs. <l> for 
each junction. 

The normalized currents i = I / Iq through the junc- 
tions are given by 



- ± j = pcik + ik + ic{^) sin(7fe) + iN.k (1) 

where /c = 1,2 labels the two Josephson junctions, /q 
is the amplitude of the Josephson current, j = J/Iq 
is the normalized circulating current in the SQUID 
loop and '±' refers to junctions 1 and 2, respectively. 
Pc = 27rIoR^C/^o is the Stewart-McCumber param- 
eter. <l>o is the flux quantum and R and C, re- 




I R=10.6mV 
J I I I 



-20 



-10 




V(mV) 



10 



20 



FIG. SI: (Color online) Current voltage {IV) characteris- 
tics of (a) the 0-design SQUID and (b) the vr-design SQUID. 
Each graph contains two measured curves and two calculated 
curves. Measurements were taken at, respectively, maximum 
(black circles) and minimum (grey circles) critical current. 
Calculated curves were obtained by numerically solving the 
SQUID Langevin equations, extended by taking the nonzero 
junction width into account, for (a) a SQUID at flux $ = 
(red line) and $ = $o/2 (magenta line), and (b) for a tt 
SQUID at flux $ = $0/2 (red line) and $ = $0 (magenta 
line). Model parameters in (a) are /3c = 0.5, /3l = 1.4, F = 
0.022, aL — 0.05 and fj — 0.128. By matching abscissa and 
ordinate of the measured and calculated curves one obtains 
/o = 8.2 /xA and IqR — 7.55 /xV. Model parameters in (b) are 
P, = 0.7, Pl = 2.2, F = 0.0147, aL = 0.05, fj = 0.12, /o = 
12.2 /i A and IqR = 10.6 /iV. The same set of parameters was 
used for all curves discussed in this paper. 



spectively, denote junction resistance and capacitance. 
7/e is the Josephson phase difference across junction k 
and the dots denote derivative with respect to time. 
The flux dependent quantity ic{^) is given by ic{^) = 



sin(7r/j^/$o)/[7r/j^/^o]- The normalized noise cur- 
rent ijy^k has a spectral power density 4r, with T = 
27rkBT/Io^O' The above equations assumed that the 
junction parameters /3c and ic{^) are the same for both 
junctions. If junction A: is a tt junction, a phase tt is to 
be added to jk- 

The two phases jk are related by 



design I 



Simulation 



72 - 71 = 27r^/^o + ttPlU + ^lO 



(2) 



where Pl = 2IqL/^q. L = Li +L2, where Li and L2 are 
the inductances of the two SQUID arms, related to the 
inductance asymmetry aL via L^ = 1/(1 ± aL)/2. 

From Eqs. (1) and (2) one obtains the normalized volt- 
age V = V/IqR^ and thus current voltage characteristics, 
critical current vs. flux etc., by taking the time average 
of (71 +72)72. A consistent set of model parameters can 
be obtained by analyzing IV characteristics at maximum 
and minimum /c, plus /c vs. H. 

The calculated IV characteristics in Fig. SI (a) are for 
/3c = 0.5, /3l = 1.4, F = 0.022, aL = 0.05 and fj = 0.128 
(the latter two numbers are actually determined from Ic 
vs. H data). In dimensioned units one finds /q = 8.2 /i A 
and IqR = 7.55 /iV, R = 0.921], L = 177 pH, C = 24 pF. 
These are reasonable numbers for our SQUIDs. In the 
/iqH = 0.1 /iT curve of Fig. SI (a) one also notices a hump 
for 6 /iV < F < 10 /iV. This is a LC resonance, which be- 
comes maximally excited when the supercurrents across 
the two junctions oscillate out-of-phase^. 

Two IV characteristics for the 7r-design SQUID are 
shown in Fig. SI (b). One first notices that the first 
Ic maximum (black circles) was obtained at a relatively 
large field, 0.07 /iT. By contrast, at jiioH — 0.01 /iT, Ic 
had a minimum (grey circles). Simulating these curves 
(red line for ii^R = 0.07 /iT and magenta line for fiQH = 
0.01 /iT) we have assumed that junction 2 is a tt junction 
and further used the parameters /3c = 0.7, /3l = 2.2, F 
= 0.0147, aL = 0.05, fj = 0.12, /q = 12.2 /iA, IqR = 
10.6 /iV, R = 0.871], C = 25pF, L = 187 pH, which are 
not very different from the reference SQUID. 



II. VOLTAGE VS. APPLIED FIELD 

Fig. S2 shows measured [(a), (c)] and calculated 
[(b), (d)] patterns V vs. H for the 0-design SQUID [(a), 
(b)] and the 7r-design SQUID [(c), (d)]. The data were ob- 
tained in the same run as the IV characteristics of Fig. SI 
and the Ic vs. H data (Fig. 2 of the main paper). In the 
measurements the current / across the 0-design SQUID 
has been increased (decreased) from in steps of 1 /iA 
up to ±28 /iA. For the 7r-design SQUID the step width 
was 0.98 /iA, with a maximum current ±39.5 /iA. Mea- 
sured data for the 0-design SQUID, cf. Fig. S2 (a), are 
smooth and well reproduced by the calculated patterns 
[Fig. S2 (b)]. Note that for voltages below |y| ^ 5 /iV 
the |y| vs. H minima are located near H = (mod- 
ulo modulation period), while for \V\ > 5/iV one finds a 




FIG. S2: (Color online) Measured V vs. H for (a) the 0- 
design SQUID and (c) the vr-design SQUID together with 
numerically calculated curves for (c) the 0-design SQUID and 
(d) the TT-design SQUID. The current / across the 0-design 
SQUID has been increased (decreased) from in steps of 1 /iA 
up to ±28 /iA. For the vr-design SQUID the step width was 
0.98 /iA, with a maximum current ±39.5 /iA. 



maximum here. This is due to the LC resonance which 
is also visible in Fig. SI (a). Further, the |U| vs. H 
maxima appear (modulo modulation period) near fioH 
= 0.995 /iT and are slightly shifted with respect to the 
origin of the ordinate. This shift, which is due to the 
small inductance asymmetry, was seen already in Ic vs. 
H. Otherwise, the device behaves as it can be expected 
for a SQUID. In the experimental data for the the tt- 
design SQUID one observes similar features, but shifted 
by a half period with respect to the 7r-design SQUID. As 
in /c vs. H, jumps appear in the second SQUID modula- 
tion period at positive H at positive V and at negative H 
at negative V. This point symmetry was already visible 
in the /c vs. H data. Apart from these jumps, V ys. H 
is reasonably well reproduced by simulations, cf. Fig. S2 
(b). 



III. CRITICAL CURRENT VS. APPLIED 
FIELD: SYMMETRY CONSIDERATIONS 

In Fig. S3 (a) we further demonstrate point symmetry 
by plotting the negative Ic as — /c vs. —/j^qH (red line) 
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FIG. S4: (Color online) Resonance current Ir vs. applied 
field for (a) the 0-design SQUID and (b) the 7r-design SQUID 
(black curves), together with the corresponding simulated 
curves (red lines). Insets show Ir vs. H at expanded cur- 
rent and field scales. 



FIG. S3: (Color online) h vs. H of the vr-design SQUID: (a) 
Test of symmetries of the Ic jumps. Plot of negative Ic as 
— /c vs. —/hqH (red line) on top of Ic vs. /j^qH for positive 
Ic (black line), demonstrating point symmetry, (b) Ic vs. H 
with an additional modulation current /mod applied across 
GBJ 1. Black line: /mod = 0; green line: /mod — — 15/iA; red 
line: /mod = 15 /xA. Inset shows first penetration fields /jqHcj 
vs. /mod- Grey line in (b) indicates offset field. 



on top of Ic vs. jiioH for the positive Ic (black line). The 
curves are basically indistinguishable. 

The point symmetry gives some hint on the proper- 
ties of the trapped flux. First, if Abrikosov vortices or 
Josephson fluxons were trapped permanently one would 
at most expect a symmetry with respect to a change 
/ ^ — /. The point symmetry implies that also the polar- 
ity of the trapped flux changes by reversing the magnetic 
field. The fact that the Ic jumps exhibit almost no hys- 
teresis further implies that almost no pinning is present. 
All this, together with the observation that jumps in Ic 
vs. H are absent for the 0-design SQUID, supports the 
suspicion that Josephson fluxons enter and leave GBJ 1. 
The change in polarity of the trapped flux with applied 
field implies that no fluxons are present at low fields, al- 
lowing to interpret our data in terms of a tt SQUID and 
thus in favor of a 6^3.2 _^2 -wave order parameter symmetry. 

The Ic jumps in Fig. 2 (b) of the main paper were 
strongly asymmetric with respect to the applied field. 
For positive current the first Ic jumps occurred at the 



penetration fields iiqHcj = +0.09 /iT and at -0.24 /iT. 
At negative current they were observable at (iqHcj = 
-0.09 /iT and at +0.24 /iT. This asymmetry was differ- 
ent in different cooldowns and seemed to depend on the 
residual field. In fact, screening currents across GBJ 1 
will change the surface barrier for fluxon entry, making 
//cj asymmetric. Fig. S3 (b) shows results from another 
cooldown where the jumps occurred almost symmetri- 
cally with respect to H (black line). In this graph we 
also show two measurements, where we applied an addi- 
tional modulation current /mod across GBJ 1, cf. Fig. 1 
in the main paper. Depending on the polarity of this cur- 
rent its effect was to linearly shift the appearance of the 
Ic jumps, i.e. Hcj^ to higher (lower) values for positive 
(negative) values of /mod- This shift did not depend on 
the sign of /. /mod exerts a Lorentz force on the fluxons. 
For /mod > fluxons appearing at // > are repelled 
from the SQUID, increasing the switching fields for both 
polarities of /. The force on antifluxons appearing at 
// < points inward, i.e., the switching field decreases. 
For /mod < the forces on fluxons and antifluxons are 
opposite. Thus, the shifts of the Ic jumps are fully com- 
patible with the notion of (anti)fluxons having entered 
GBJ 1. 



IV. LC RESONANCES 

We finally address in more detail the resonance feature 
which appeared as a hump in the IV characteristics at 
jj^qH = 0.1 /iT for the 0-design SQUID and at jj^qH = 
0.01 /iT for the 7r-design SQUID. In V vs. H the reso- 
nance appeared as a phase shift of a half period in the 
SQUID modulation. These effects are due to resonantly 
excited circulation currents across the SQUID loop and 
occur when the Josephson currents across the two junc- 
tions acquire an out-of phase component. Maximum exci- 
tation occurs when the Josephson currents oscillate maxi- 
mally out-of-phase. For a SQUID one thus expects the 
strongest effect for a half flux quantum applied to the 



junction, for a tt SQUID it should be strongest at integer 
multiples of <I>o^. To investigate the LC resonance vs. 
applied field H we have determined the "resonance cur- 
rent" Ir across the SQUIDs by using a voltage criterion 
Vc which corresponds to the center of the hump feature 
in the IV characteristics (7.5 /iV for the 0-design SQUID 
and 7/iV for the 7r-design SQUID). 

Fig. S4 (a) shows Ir vs. H for the 0-design SQUID 
(black line) together with a simulated curve (red line). 
Apart from a slight offset due to residual fields the mini- 
mum in Ir is at zero field both for negative and positive 
values of Ir. By contrast, for the 7r-design SQUID one 
finds (except for an offset) a maximum of Ir near zero 
field, again indicative of a tt SQUID. 
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